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Abstract 

We report on the production of J/ip mesons in 530 and 800 GeV/c pBe collisions in 
the Feynman-x range 0.0<a;i;'<0.6. The J/ip mesons were detected via decays into fi'^fi". 
Differential distributions for J/ip production have been measured as functions of xp, Pti 
and cosine of the Gottfried- Jackson decay angle. These distributions are compared with 
results on J /ip production obtained in 515 GeV/c n^Be collisions, measured by the same 
experiment, as well as with results from other experiments using incident protons. 

PACS numbers 13.85.Ni, 12.38.Qk, 25.40.Ve 
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I. INTRODUCTION 

Studies of charmonium production in hadron collisions provide important information 
on both perturbative and non-perturbative QCD. Recent advances in the understanding 
of quarkonium production have been stimulated by the unexpectedly large cross sections 
observed for direct J /ip and ip{2S) production at large px measured at the Fermilab Teva- 
tron ||T]]. These observations have led to renewed interest in this field and to an improved 
theoretical understanding of J/ip production mechanisms. In particular, a new model has 
been developed to explain the CDF results: the color-octet mechanism [Q,^. In the original 
color-singlet model @], charmonium mesons retain the quantum numbers of the cc pairs 
produced in the hard scatters, and thus each J^*-^ state can only be directly produced via 
the corresponding hard scattering color-singlet sub-processes. The color-octet mechanism 
extends the color-singlet approach by taking into account the production of cc pairs in 
a color-octet configuration accompanied by a gluon. The color-octet state evolves into a 
color-singlet state via emission of a soft gluon. 

Another recent analysis of J/ip photo- and hadro-production data resurrects the 
twenty-year-old color evaporation approach. This approach differs from Ref. 0,^ in the way 
the cc pair exchanges color with the underlying event. It assumes multiple exchanges of soft 
gluons, whereas such interactions are suppressed in the color octet model by powers of v, 
the relative velocity of the heavy quarks within the cc system. These multiple soft gluon 
exchanges destroy the initial polarization of the heavy quark pair. Therefore, polarization 
measurements may provide the best tool to distinguish between these phenomenological 
models. 

In addition to the Tevatron collider data, a number of new results have become available 
at fixed-target energies fl^TO| . These measurements provide further tests for the color octet 



and color evaporation models 0. 

In this paper, we report on differential distributions for J/ip production in 530 GeV/c 
and 800 GeV/c pBe collisions, and compare these distributions with our published results 



based upon studies of 515 GeV/c t:~ Be collisions 0. The pion and proton data samples were 
collected during the 1990 and 1991 Tevatron fixed-target runs, respectively, using essentially 
the same apparatus. In addition, during the 1991 fixed-target run we recorded a lower 
statistics sample of data using a 515 GeV/c negative pion beam to verify the consistency of 
data samples from the two runs. 

We also compare the distributions from the 800 GeV/c incident proton data with the 
recently published results from Fermilab experiments E789 PJ^ and E771 [l^ obtained at 
the same energy. The xp range covered by this analysis (0.0<xf<0.6) complements the xp 



ranges studied in Ref. g] (0.3<a;F<0.95), in Ref. || (-0.035<2;i7<0.135), and in Ref. [ig 
(-0.05<xf<0.25). 

II. APPARATUS 



The experiment was performed in the Fermilab Meson West beam line using a large- 
aperture, open-geometry spectrometer with the capability of triggering on the presence of 



high-mass muon pairs P,p]l|. The layout of the experiment is shown in Fig. 1. A brief 
description of the essential elements of the Meson West apparatus as implemented during 
the 1991 fixed-target run is presented below. 

The beam line included a 42-m-long differential Cerenkov counter capable of tagging 
incident beam particles. The unseparated positive secondary beam at 530 GeV/c was com- 
posed primarily of protons with a small admixture (<3%) of pions [l^. The 800 GeV/c 
beam was a coUimated primary proton beam extracted from the Tevatron. The beam in- 
tensity varied during the run, and at its maximum was 2x10® protons per 23-second-long 
spill (at 57 second intervals). 

A hadron shield consisting of 5 m of iron was located upstream of the target to reduce 
background from off-axis hadrons and soft muons. Two scintillator veto walls were placed 
upstream of the hadron shield, and two others were located downstream of the shield. These 
veto walls were used to identify high energy muons that penetrated the iron. An on-line veto 



in the dimuon trigger rejected events containing coincidences between signals from either of 
the upstream veto walls and either of the downstream veto walls. 

The incident beam was defined via a three plane (X, Y, and U) beam hodoscope with 
millimeter-scale resolution, combined with a beam halo (BH) system of four scintillation 
counters located upstream of the target. A beam particle was identified by at least a two- 
fold coincidence from the beam hodoscope in the absence of any signal from the BH system. 
The BH counters formed a rectangle with a 1 cm diameter hole centered on the beam. 
Signals from the beam hodoscope were also employed to reject events containing more than 
one spatially isolated incident beam particle. 

The target consisted of two copper pieces, each 0.8 mm thick, followed by a 15-cm- 
long liquid hydrogen target with 0.25 cm and 0.28 cm thick beryllium windows upstream 
and downstream of the hydrogen target, respectively. A 2.54 cm thick beryllium piece was 
located immediately downstream of the second window. The total thickness of the target 
material corresponded to ^11% of a proton interaction length. 

Three X-Y modules of silicon-strip detectors (SSDs) were located upstream of the target 
to measure the trajectories of incident beam particles. Each module consisted of two single- 
sided planes. Five more X-Y SSD modules were located downstream of the target. These 
modules were used to measure the trajectories of charged particles produced in the target 
and to reconstruct primary and secondary vertices. The SSDs in the first downstream 
module had a 25-/im-pitch inner region and a 50-/xm-pitch outer region; the SSDs in the 
other modules had 50 /im pitch throughout. 

A dipole magnet, producing a px impulse of 0.45 GeV/c, was located downstream of 
the SSDs. Four proportional wire chamber (PWC) modules were located downstream of 
the dipole magnet. Each module contained four wire planes to provide X, Y, U, and V 
measurements. (The U and V coordinates were a pair of orthogonal coordinates rotated 
by 37° relative to the X and Y coordinates.) The wires in each plane had 2.54 mm pitch. 
An area of from 6 to 26 cm^ in the center of each plane had reduced sensitivity to provide 
protection against beam particles. Two stations of straw tube drift chambers, each with four 



planes in the X and four planes in the Y directions, were located adjacent to the outside 
PWC modules. In the analysis reported here, momentum measurements for charged tracks 
employed all applicable tracking devices: the SSDs, the dipole magnet, the PWCs, and the 
straw tubes, yielding a momentum resolution of Ap/p = 0.02% p. In addition, the muon 
track candidates were constrained to come from the primary vertex. (The straw tube data 
were not used in the results reported in Ref. 0.) 

Interactions were detected by pairs of scintillation counters located upstream and down- 
stream of the dipole magnet. An interaction was defined by signals from at least two of these 
four counters in coincidence with a valid beam particle. The interaction rate was generally 
less than 1 MHz. 

A liquid-argon calorimeter (LAC) was located downstream of the magnetic spectrometer. 
The LAC contained both electromagnetic and hadronic sections. The outer radius of the 
electromagnetic section was 165 cm; the inner radius was 20 cm. A helium-gas-filled beam 
pipe was inserted along the axis of the LAC to minimize interactions of beam particles. The 
total LAC material corresponded to more than 10 absorption lengths. 

An iron and scintillator forward calorimeter was located downstream of the LAC to detect 
particles passing through the beam pipe. The forward calorimeter had an outer diameter 
of 1 m with a 3 cm diameter axial beam hole, and represented 10 absorption lengths of 
material. 

The muon detector was located 20 m downstream of the target and extended for 16 m. 
The detector contained, in sequence, two muon PWC stations (/xA, /iB) with four planes 
each (X, Y, U, and V), a beam dump consisting of tungsten and steel imbedded in concrete, 
an iron toroid magnet producing an average pt impulse of 1.3 GeV/c, and four more PWCs, 
(yul — yu4), each with three planes (X, U, and V). In these PWCs, the (7 and V coordinates 
were at angles of 45° above and below the horizontal, respectively. Iron, lead, and concrete 
shielding was interspersed between chambers fil through /i4. Two muon hodoscope planes, 
HI and H2 (each with sixteen petal-shaped scintillation counters), were also located in 
this region. The outer radius of the muon chambers and muon hodoscopes was 135 cm. 



The acceptance of the muon spectrometer was hmited by tapered axial holes, ranging from 
13 cm to 20 cm in radius, through the toroid magnet, the muon chambers, and the muon 
hodoscopes. The hole in the toroid was filled with lead absorber. 

The muon detector elements HI, H2, and /il through /i4 were shielded from hadrons 
by the material in the LAC, the forward calorimeter, the beam dump, the toroid, and 
the concrete shielding. Muons produced at the target with energies larger than ^15 GeV 
penetrated this material, and all particles reaching the muon hodoscopes were assumed to 
be muons. Two or more hits in each of the muon hodoscopes were required as part of 
the dimuon pretrigger; the average hit multiplicity in dimuon triggered events was 2.3 per 
hodoscope plane. In addition to the muon hodoscope requirement, the upstream veto wall 
requirements completed the pretrigger conditions. The pretrigger rate was 1.2 xlO^^ and 
2 xlO~"^ per live interaction, for 530 GeV/c and 800 GeV/c incident protons, respectively. 
No radial dependence of the pretrigger counter efficiencies was observed. 

Events satisfying the dimuon pretrigger were evaluated by the dimuon trigger processor 
(DMTP), which reconstructed space points in PWCs /il and /i4, formed muon tracks (re- 
quiring an additional hit in either /i2 or fi3 along the projected track trajectory), measured 
muon momenta from the estimated bend in the toroid (assuming that the tracks originated 
in the target), and calculated the dimuon invariant mass. Trajectories of muon tracks re- 
constructed in the downstream muon PWCs were projected to the center of the target. The 
processor accepted only muon tracks with at least three hits within roads around the pro- 
jected trajectory in both the fiA and /iB chambers. A dimuon trigger resulted if any of the 
calculated dimuon masses was above a preset threshold. A mass threshold of 1.0 GeV/c^ 
resulted in a trigger rate of 1 xlO^^ to 2 x 10^^ per live interaction. The average DMTP 
processing time was 10 fis per pretrigger, which included 5 fis to decode the muon chamber 
data. The J/ip mass resolution of the DMTP was 550 MeV/c^. 

The toroid polarity was reversed several times during the run. In addition, part of the 
data was recorded without the /lA and /iB hit requirement in the trigger. 



III. DATA 

Table I lists the numbers of dimuon triggers and approximate integrated luminosities per 
nucleon on Be recorded during the 1991 fixed-target run for the different incident beams 
employed. For each dimuon-triggered event, the reconstructed muon tracks were linked 
through the entire detector. Only events with at least two fully-linked muons were processed 
further. The remaining track segments in the SSDs and upstream PWCs were used in the 
reconstruction of other tracks and event vertices. The distribution of the reconstructed 
primary vertices along the nominal beam direction (Z), for events containing reconstructed 
dimuons in the mass range 2.8 GeV/c^ to 3.4 GeV/c^, is shown in Fig. 2. 

Dimuons contributing to this analysis came from events with primary vertices in the 
beryllium targets, and had dimuon Feynman-a; [xp = 2^2/1/5) in the range 0.0<a;i;'<0.6 
(0.0< xf<0.8 for 7r~ data). Data from other targets are used to investigate nuclear effects. 
Figure 3 shows the reconstructed opposite-sign dimuon invariant mass distributions in the 
J/ip mass region for the various data samples. The full width at half maximum (FWHM) 
for the J/ip signals varied from 120 to 135 MeV/c^. Fits were performed on the dimuon 
mass spectra in the mass range 1.5 to 5.0 GeV/c^. These fits employed the sum of two 
exponentials for the continuum background, and resolution functions for the J ftp and ip{2S) 
resonances determined via Monte Carlo simulations. The J/ip masses resulting from the fits 
are within 1 MeV/c^ of the Particle Data Group (PDG) value [^. The ip{2S) mass was 
fixed at the PDG value. The numbers of J /ips and ip{2S)s obtained from these fits are listed 
in Table I. The systematic uncertainties in the number of J /ip and 'ip{2S) combinations take 
into account variations in the shapes assumed for the signal and background distributions. 

To evaluate reconstruction efficiencies and acceptances, 5 x 10^ Monte Carlo events con- 
taining J/ips were generated. The xp and p"^ distributions of the generated Monte Carlo 
J/ips were based on our previous measurements 0,0. We assumed that the J/ips decay 



isotropically (consistent with observations to be described in the next section). These Monte 
Carlo events also contained additional charged tracks with an average multiplicity consistent 
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with the data. The dimuons and associated particles were propagated through a GEANT 
simulation of the Meson West spectrometer, which incorporated measured spectrometer 
chamber efficiencies and instrumental noise determined from our data. The dimuon events 
were then reconstructed using the same tracking programs used for reconstruction of the 
data. We have evaluated the J/ip reconstruction efficiency and geometrical acceptance as a 
function of three variables: (i) xp-, the J /ip Feynman-x; (ii) p^, the square of the J /ip trans- 
verse momentum; and (iii) cos9, the cosine of the Gottfried- Jackson decay angle between the 
/i"*" and the beam axis in the J /ip rest frame. The products of acceptance and reconstruction 
efficiency, aji^^^+^- ■ eji^^^+^- , were evaluated as one-dimensional distributions averaged 
over the other variable using an iterative technique (see Fig. 4). For a given iteration, 
the input Monte Carlo xp and p\ distributions were weighted to match those of the data 
corrected by the acceptance and the reconstruction efficiencies determined in the previous 
iteration. The iterations were halted once stability was achieved. In preparing Ref. we 
demonstrated that this method of evaluating acceptance and efficiency yielded cross section 
results consistent with a more sophisticated method based on a two-dimensional acceptance 
surface over the xp and p\ plane. The resulting xp and p\ acceptances and efficiencies 
are almost independent of the toroid polarity. But the product of acceptance and efficiency 
as a function of cosO, integrated over the xp and p\ spectra, shown in Fig. 4, exhibits a 
significant dependence on the toroid polarity. This is a result of the difference in acceptance 
as muons are either "bent in" or "bent out" by the toroid magnet in the muon detector. 

IV. J/i) DIFFERENTIAL DISTRIBUTIONS 

J/ip candidates consisted of opposite sign dimuons with invariant mass between 
2.8 GeV/c^ and 3.4 GeV/c^ originating from the primary vertex. Fits to the dimuon mass 
spectra outside the resonance band, similar to the ones described in Section III and shown 
in Fig. 3, for different regions of xp, p"^, and cos6 indicate significant variations in the back- 
ground contribution as a function of kinematic variables (for the 800 GeV/c proton data). 
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Therefore, in contrast to our previous analysis P], we have not carried out a two-constraint 
kinematic fit for each J/ip candidate to improve the muon momentum resolution. Instead, 
we performed a bin-by-bin background subtraction to measure the J/ip differential distri- 
butions. The cosO distributions were measured separately for the positive and negative 
toroid polarities, and averaged as independent measurements. The effects of experimental 
resolution were taken into account in the acceptance corrections. 

The xp distributions for the incident proton data at 530 GeV/c include corrections for 
the 2.8% incident vr"^ contamination [|l^]. However, the Cerenkov counter information was 



not available for the entire data set. In addition, the Cerenkov tag for protons that results 
in 99.6% protons among the selected tagged particles also results in ^ 25% reduction in the 
available luminosity. Therefore, rather than use the beam Cerenkov counter information on 
an event-by-event basis, we applied iterative statistical corrections to individual xp bins. 
The corrections were based on the measured J /ip integrated cross sections per nucleon times 
J 1'^-^ fJ-^fJ'^ branching ratio for xp>0 of 9.2 ± 2.0 nb |I^ and 12.9 ± 1.6 nb [^ for incident 



protons and pions, respectively. The shape of the xp distribution for incident pions was 
taken from Ref. ^. The size of the correction varied from 1% at Xi;'=0.0 to 22% at xp= 
0.55. For a given xp bin, the systematic uncertainties due to this correction are much 
smaller than the corresponding statistical uncertainties. The effects of beam contamination 
are negligible for the p^, and cos6 distributions, which are similar for incident pions and 
protons. 

The J/ip differential distributions are shown in Figs. 5 through 7, and tabulated in Ta- 
bles II and III. Data points were corrected for the geometrical acceptances and reconstruction 
efficiencies discussed in the previous section. The quoted uncertainties represent statistical 
and systematic uncertainties added in quadrature. The systematic uncertainties include 
contributions due to uncertainties in the background estimates, as well as pion contamina- 
tion for the 530 GeV/c incident proton data sample. Additional systematic uncertainties, 
due to the trigger and off-line reconstruction acceptance and efficiency calculations, were 
estimated from the variation of our results for data subsamples taken under different run- 
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ning conditions. These contributions are included only in the systematic uncertainties of 
the parameters of empirical fits employed to describe the data. Results of these empirical 
fits are summarized in Tables IV through VI. 
We used the parameterization: 

da/dxp oc (1 — xfY, (4.1) 

to describe the J/tp xp distribution for the proton-induced collisions [see Figs. 5(a) and 
6(a)]. The fit value of parameter c increases from c = 5.89 ± 0.14 (stat) ± 0.11 (syst) 
at 530 GeV/c to c = 6.18 ± 0.12 {stat) ± 0.11 (syst) at 800 GeV/c. Our resuhs are 
consistent with the general trend of increasing c as the center-of-mass energy increases (see 
Fig. 8). However, the values of parameter c for the incident beam momentum of 800 GeV/c, 



quoted by other experiments, are 6.38±0.24 for a silicon target [|10| and 4.91±0.18 for a gold 
target 0. (The silicon target data is restricted to the xp range -0.05<xf<0.25, while the 
gold target data is restricted to -0.035<Xi7<0.135). In Fig. 9(a), we directly compare the xp 
distributions from the three experiments, normalized to the integral over the 0<Xi7'<0.135 
range. Our results agree very well with those of Ref. [|r^, whereas the data of Ref. 
exhibit a harder xp dependence. Quantatively, when a fit of the xp distribution for our 
800 GeV/c data is restricted to the 0<Xi7'<0.15 range, that fit yields c = 5.94 ± 0.43 with 
a x^ psr degree of freedom, x^/^df = 6.6/7. Previous studies of the nuclear dependence 
of J/ip production have exhibited little nuclear dependence at small xp [jl5[. Comparisons 
of J/ip Xp distributions from our data using H2,Be, and Cu targets, presented in Fig. 10, 
do not show any significant dependence on the target nucleus. A difference of one unit 
in the value of the parameter c implies a relative yield change between xp= 0.0 and xp= 
0.25 by a factor of 0.75. The data show that the ratios of the normalized distributions for 
different nuclear targets are consistent with unity (within 25%) over the entire measured 
Xp range. Our observation of weak xp dependence for nuclear effects (within our large 
statistical uncertainties) is consistent with results for pCu and pBe collisions (xi?>0.3) 
reported in Ref. |^. The xp distributions from Ref. [§ are compared to our measurements 

11 



in Fig. 9(b). Reasonable agreement is exhibited in the overlapping xp region. We suspect 
that the discrepancy observed in the small xp region (discussed above) is due either to a 
large statistical fluctuation or neglected experimental effects in some of the data sets. 



We have also used the parameterization proposed in Ref. [|I^ to describe da/dxp '■ 

da/dxp oc (1 — Xi)'^(l — X2)'^/{xi + X2), (4.2) 

where 



xi,2 = 0.5{^Jxf^ + 4:T±xf), (4.3) 

and T = Mj,^/s. Here, xi represents the beam parton fractional momentum and X2 the 
target parton fractional momentum. The values for the parameter k, obtained from the fit 
are 4.95 ± 0.17 (stat) ± 0.12 (syst) at 530 GeV/c and 4.95 ± 0.16 (stat) ± 0.12 (syst) 
at 800 GeV/c. The corresponding x^ /ndf are 24/13 and 30/13. A similar fit to the vr" 
induced data at 515 GeV/c 0, with the (1 — a;2)'" factor replaced by {1 — X2)'^^ , gave 
K=1.69 ± 0.04 with 'X^ /ndf = 50/33. According to Ref. |T^, if gluon fusion is the dominant 
mechanism responsible for J/tp production via incident protons, then the expected values of 
the parameter k are 5 for J/ips produced predominantly through x states decays, and 6 for 
direct J/ip production. For incident pions, the gluon fusion mechanism results in k values 
of 3 and 4 for the indirect and direct J/tp production, respectively. The qq annihilation 
mechanism results in values of k smaller by one unit. Our data indicate a larger difference 
between J/ip production by incident protons and pions than predicted by simple counting 
rules. 

Gupta and Sridhar made a detailed comparison of the color-octet model predictions 
with measured xp distributions from several fixed-target experiments, including our previous 
measurements described in Refs. §,0]. AH values of the required non-perturbative matrix 
elements were derived from other experiments, primarily from charmonium photoproduction 
data and from J/ip production data from the Tevatron. In general, the parameter-free 
predictions of the model for XF-distributions agree with the proton and pion induced data 
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both in magnitude and shape. However, uncertainties in the model predictions due to scale 
uncertainties, the various parton density sets, and non-perturbative matrix elements are still 
larger than the experimental uncertainties. 



The p^ spectra (Figs. 5b and 6b) were fit using the empirical form |]T7 



da/rfp^oc(l+p^//?^)-^ (4.4) 

with parameters /3 = [2.66 ± 0.04 {stat) ± 0.03 {syst)] GeV/c and [2.81 ± 0.03 {stat) ± 
0.03 {syst)] GeV/c for the 530 GeV/c and 800 GeV/c data, respectively. The corresponding, 
XiT'-integrated average J /ip transverse momenta, (pt), are: [1.15 ± 0.02 {stat)] GeV/c and 
[1.22 ± 0.01 {stat)] GeV/c. As shown in Fig. 11, the value of {px) grows slowly with center- 
of-mass energy. The comparison of J /ip pt distributions for the H2, Be, and Cu targets, 
presented in Fig. 12, does not show any significant dependence on the target nucleus. 

Fits to the cos9 distribution in the range \cos9\<0.8 [Figs. 5(c) and 6(c)] using the 
functional form 

da/d{cose) oc{l + Xcos'^e), (4.5) 

yield A = 0.01 ± 0.12 {stat) ± 0.09 {syst) and A = -0.11 ± 0.12 {stat) ± 0.09 {syst) at 
530 GeV/c and 800 GeV/c, respectively. Figures 5(c) and 6(c) show the data for each toroid 
magnet polarity separately, and the results are similar. As illustrated in Fig. 4, the accep- 
tance and efficiency as a function of cos6 are sensitive to the toroid polarity. Nevertheless, 
all data sets are consistent with unpolarized J/ip production. The E771 Collaboration re- 



ported A = -0.09±0.12 |T^ in pSi interactions at 800 GeV/c. Our measurement of A for 



J/ip production from the high-statistic incident vr data sample gave A = -0.01±0.08 [Q. 



These results confirm some of the earlier observations which reported small values of A [jl8 
and are consistent with expectations from the color evaporation model [0. These results 
for A are in disagreement with expectations from the color octet model ||^ , which predicts a 
sizeable J/ip transverse polarization, with mild energy dependence (e.g. 0.31 < A < 0.63 at 
v^= 15.3 GeV). 
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The differential distributions for J /ip production via incident 515 GeV /c vr" beams are 
shown in Fig. 7. The new, lower statistics results, obtained during the same run as the 
proton data, are in a good agreement with our published high-statistics results . The only 
substantial differences between the J /ip differential distributions from the pion and proton 
induced data samples are in the xp distributions. These differences are attributed to the 
harder distribution of gluons within the incident pion as compared to an incident proton. 

V. SUMMARY 

We have studied the production of J /if) mesons in the Feynman-x range 0.0<a;i;'<0.6 
in 530 GeV/c and 800 GeV/c pBe collisions. We have parametrized the J /ip differential 
distributions and measured their dependence on center-of-mass energy. The xp distributions 
of J /if) production become more central as ^/s increases. For incident 800 GeV/c protons, 
the Xp distribution observed in our data {Be target) is similar to the distribution reported 



by E771 {Si target) ||T^, and differs from the distribution reported by E789 {Au target) 
0. Based upon studies of data from H2,Be, and Cu targets, we are inclined to believe 
that this difference is not likely to be explained by nuclear effects. The xp distribution 
for J/ip production in tt~ Be interactions is less central than the corresponding distribution 
for pBe interactions at a similar beam energy. The data are consistent with unpolarized 
J /ip production as anticipated in the color evaporation model, but are in disagreement with 
expectations from the color octet model. 

We thank the staffs of all participating institutions, especially those of Fermilab and the 
Institute for High Energy Physics (IHEP) at Protvino. This work was supported by the 
U. S. Department of Energy, the National Science Foundation, and the Russian Ministries 
of Science and Atomic Energy. 
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Table I. Data sets and number of J /ip and ip{2S) combinations from interactions on the 

Be targets during the 1991 fixed-target run. 



Beam 


Momentum 


Integrated 
Luminosity 


Dimuon 

triggers 


Nj/^ 


N^i2S) 




(GeV/c) 


(pb-i) 


103 






TT^ 


515 


1.4 


902 


806±32 ± 13 


18 ±10 ±6 


P 


530 


6.4 


2434 


3607±65 ± 28 


100 ±18 ±15 


P 


800 


7.3 


6400 


6540±90 ± 43 


214 ±27 ±20 
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Table II. Differential distributions for J /ip production in 530 GeY /c pBe interactions. The 

quoted uncertainties represent statistical and systematic uncertainties added in 
quadrature. The systematic uncertainties include contributions due to uncertainties in the 

background estimates and pion contamination. 



Xf 




{l/a)d(r/dxF 


Pt 

(GeV/c)2 


{l/(x)d(T/dpl 

l/(GeV/c)2 


cos9 




{l/a)da/dcos9 


0.00- 


-0.05 


5.7 ±1.6 


0.00-0.25 0.659 ±0.028 


-0.8- 


-0.6 


0.532 ±0.038 


0.05- 


-0.10 


4.56 ±0.34 


0.25-0.50 0.556 ±0.026 


-0.6- 


-0.5 


0.490 ±0.027 


0.100- 


-0.125 


3.60 ±0.25 


0.50-0.75 0.401 ±0.023 


-0.5- 


-0.4 


0.523 ±0.024 


0.125- 


-0.150 


2.90 ±0.19 


0.75- 1.00 0.311 ±0.021 


-0.4- 


-0.3 


0.463 ± 0.020 


0.150- 


-0.175 


2.43 ±0.15 


1.00- 1.25 0.305 ±0.021 


-0.3- 


-0.2 


0.507 ±0.018 


0.175- 


- 0.200 


1.92 ±0.12 


1.25- 1.50 0.250 ±0.019 


-0.2- 


-0.1 


0.522 ±0.020 


0.200 - 


- 0.225 


1.93 ±0.11 


1.5-2.0 


0.187 ±0.012 


-0.1- 


0.0 


0.502 ±0.021 


0.225 - 


-0.250 1.183 ±0.086 


2.0 - 2.5 


0.125 ±0.010 


0.0- 


-0.1 


0.539 ±0.022 


0.250 


-0.275 1.185 ±0.083 


2.5 - 3.0 


0.105 ±0.009 


0.1 - 


-0.2 


0.485 ± 0.024 


0.275 - 


-0.300 0.860 ±0.071 


3.0 - 4.0 


0.064 ±0.005 


0.2- 


-0.3 


0.480 ±0.027 


0.30- 


0.35 


0.721 ±0.045 


4.0 - 5.0 


0.053 ±0.005 


0.3- 


-0.4 


0.523 ±0.032 


0.35- 


-0.40 


0.452 ±0.036 


5.0 - 7.0 


0.018 ±0.002 


0.4- 


-0.5 


0.423 ± 0.038 


0.40- 


-0.45 


0.234 ±0.027 


7.0 - 10.0 


0.006 ±0.001 


0.5- 


-0.6 


0.543 ± 0.040 


0.45- 


-0.50 


0.167 ±0.023 






0.6- 


-0.8 


0.96 ±0.214 


0.50- 


-0.60 


0.095 ±0.014 
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Table III. Differential distributions for J/ip production in 800 GeV /c pBe interactions. 
The quoted uncertainties represent statistical uncertainties and systematic uncertainties in 

the background estimates added in quadrature. 



Xf 




{l/a)da/dxF 


Pt 

(GeV/c)2 


{l/a)da/dp\ 
l/(GeV/c)2 


cosd 




{l/a)da/dcos9 


0.00- 


-0.05 


6.14 ±0.33 


0.00-0.25 0.648 ±0.021 


-0.8- 


-0.6 


0.444 ± 0.038 


0.05- 


-0.08 


4.79 ±0.22 


0.25-0.50 0.481 ±0.019 


-0.6- 


-0.5 


0.488 ± 0.030 


0.075 - 


-0.100 


4.07±0.17 


0.50-0.75 0.369 ±0.017 


-0.5- 


-0.4 


0.503 ±0.026 


0.100- 


-0.125 


3.56 ±0.14 


0.75- 1.00 0.337 ±0.016 


-0.4- 


-0.3 


0.453 ± 0.024 


0.125- 


-0.150 


2.87±0.12 


1.00- 1.25 0.289 ±0.015 


-0.3- 


-0.2 


0.476 ± 0.022 


0.150- 


-0.175 


2.29 ±0.10 


1.25- 1.50 0.226 ±0.013 


-0.2- 


-0.1 


0.500 ±0.022 


0.175- 


-0.200 1.962 ±0.094 


1.5-2.0 


0.180 ±0.009 


-0.1- 


0.0 


0.499 ± 0.023 


0.200 - 


-0.225 1.710 ±0.085 


2.0 - 2.5 


0.144 ±0.008 


0.0- 


-0.1 


0.535 ±0.023 


0.225 - 


-0.250 1.579 ±0.076 


2.5 - 3.0 


0.111 ±0.007 


0.1- 


-0.2 


0.501 ±0.024 


0.250 


-0.275 1.062 ±0.071 


3.0 - 4.0 


0.069 ±0.004 


0.2- 


-0.3 


0.493 ±0.027 


0.275 - 


-0.30 


0.837 ±0.067 


4.0 - 5.0 


0.054 ±0.003 


0.3- 


-0.4 


0.562 ±0.033 


0.30- 


0.35 


0.594 ±0.043 


5.0 - 6.0 


0.032 ±0.003 


0.4- 


-0.5 


0.468 ± 0.039 


0.35- 


-0.40 


0.341 ±0.037 


6.0 - 8.0 


0.013 ±0.001 


0.5- 


-0.6 


0.553 ±0.054 


0.40- 


-0.45 


0.305 ±0.036 


8.0 - 10.0 


0.007 ±0.001 


0.6- 


-0.8 


0.74 ±0.116 


0.45- 


0.55 


0.128 ±0.026 
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Table IV. Results of fits to the da/dxp distributions for J/ip production. 



Beam 


Momentum (GeV/c) 


c 


xVndf 


K 


xVndf 


TT^ 


515 [6] 






1.69 ±0.04 


50/33 


P 


530 


5.89 ±0.14 ±0.11 


21/13 


4.95±0.17±0.12 


24/13 


P 


800 


6.18 ±0.12 ±0.11 


21/13 


4.95±0.16±0.12 


30/13 



Table V. Results of fits to the da/dp'^ distributions for J/ip production. 



Beam 


Momentum (GeV/c) 


P 


x'/ndf 


(Pt) (GeV/c) 


71- [6] 


515 






1.17 ±0.02 


P 


530 


2.66 ±0.037 ±0.03 


19/11 


1.15 ±0.02 


P 


800 


2.81 ±0.029 ±0.03 


35/12 


1.22 ±0.01 



Table VI. Results of fits to the da/d{cos6) distributions for J /ip production. 



Beam 


Momentum (GeV/c) 


A 


xVndf 


TT^ [6] 


515 


-0.01 ±0.08 


33/34 


P 


530 


0.008 ±0.12 ±0.09 


20/12 


P 


800 


-0.11 ±0.12 ±0.09 


19/12 
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FIGURE CAPTIONS. 

Figure 1. Plan view of the Fermilab Meson West spectrometer as configured for the 
1991 fixed-target run. 

Figure 2. Primary vertex Z-coordinate distribution for events containing dimuons with 
reconstructed masses in the J /ip mass range. Unmarked spikes in the distribution 
upstream of the Cu target and downstream of the Be target are due to interactions 
in the sihcon-strip detectors. 

Figure 3. The invariant mass distributions of fi^fi~ pairs for: (a) 530 GeV/c inci- 
dent protons, (b) 800 GeV/c incident protons, and (c) 515 GeV/c incident vr". The 
sohd curve in each plot is a fit to the data; the dashed curve shows the background 
contribution. 

Figure 4. Product of geometrical acceptance and reconstruction efficiency for detection 
of J/ip — ^ fi~^fi~ as a function of: (a,b) xp, (c,d) p^, and (e,f) cos6 for 530 GeV/c (left 
column) and 800 GeV/c (right column) incident protons. The solid and dashed curves 
represent the results for the positive and negative toroid polarity, respectively. 

Figure 5. Differential distributions for J/ip production as functions of: (a) xp, (b) p'^ 
((GeV/c)^), and (c) cos6 for pBe interactions at 530 GeV/c. The cos6 distributions 
are shown separately for each toroid magnet polarity. The integrals of the distributions 
in (a) and (b) are normalized to unity. The integrals of the cos6 distributions within 
the -0.6< cos9 < 0.6 range are normalized to 0.6. The solid and dashed curves in (a) 
represent empirical fits to the data using Eqs. (4.1) and (4.2), respectively. The solid 
curves in (b) and (c) represent fits using the functions shown in Eqs. (4.4) and (4.5), 
respectively. 

21 



Figure 6. Differential distributions for J /ip production as functions of: (a) xp^ (b) p^ 
((GeV/c)^), and (c) cosO for pBe interactions at 800 GeV/c. Tlie cosO distributions 
are sliown separately for each toroid magnet polarity. The integrals of the distributions 
in (a) and (b) are normalized to unity. The integrals of the cosO distributions within 
the -0.6< cos9 <0.6 range are normalized to 0.6. The solid and dashed curves in (a) 
represent empirical fits to the data using Eqs. (4.1) and (4.2), respectively. The solid 
curves in (b) and (c) represent fits using the functions shown in Eqs. (4.4) and (4.5), 
respectively. 

Figure 7. Differential distributions for J/if) production as functions of: (a) xp-, (b) 
p\ ((GeV/c)^), and (c) cos9 for n^Be interactions at 515 GeV/c. The integrals of 
the distributions in (a) and (b) are normalized to unity. The integral of the cos6 
distribution within the -0.6< cos6 < 0.6 range is normalized to 0.6. The data of 
Ref. are shown for comparison. 

Figure 8. Center-of-mass energy dependence of the exponent c in the parametrization 
of the xp distributions for J/ip production in nucleon-nucleus interactions. The xp 
distributions for Refs. [20,21,22] were converted into the da/dxp form from the quoted 
Lorentz invariant cross sections Eda/dxp using the average value of the J/ip transverse 
momentum at a given center-of-mass energy. 

Figure 9. Differential distribution for J/ip production in 800 GeV/c pBe interactions 
as a function of xp compared with: (a) the E771 (pSi) [10] and E789 (pAu) [9] results 
(each data sample is normalized by an integral of the distribution over the 0.0 < xp < 
0.135 range); (b) the E789 (pBe) and {pCu) [8] results (each data sample is normalized 
by an integral of the distribution over the 0.3 < xp < 0.55 range). 

Figure 10. Ratios of differential distributions as functions of xp for J/ip production 
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in: (a) pp and pBe interactions, and (b) pCu and pBe interactions. The open and full 
circles represent the 530 GeV/c and 800 GeV/c data, respectively. The integrals of 
the input distributions are normalized to unity. 

Figure 11. Center-of-mass energy dependence of the average J/ip transverse momen- 
tum for nucleon-nucleus interactions. 

Figure 12. Ratios of differential distributions as functions of p"^ for J/ip production 
in: (a) pp and pBe interactions, and (b) pCu and pBe interactions. The open and full 
circles represent the 530 GeV/c and 800 GeV/c data, respectively. The integrals of 
the input distributions are normalized to unity. 
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FIG. 2. Primary vertex Z-coordinate distribution for events containing dimuons with recon- 
structed masses in the J/^ mass range. Unmarked spikes in the distribution upstream of the Cu 
target and downstream of the Be target are due to interactions in the silicon-strip detectors. 
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FIG. 4. Product of geometrical acceptance and reconstruction efficiency for detection of 
J/ip -^ fJ'~^fJ'~ as a function of: (a,b) xp, (c,d) p'^, and (e,f) cos9 for 530 GeV/c (left column) and 
800 GeV/c (right column) incident protons. The solid and dashed curves represent the results for 
the positive and negative toroid polarity, respectively. 
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FIG. 5. Differential distributions for J/tp production as functions of: (a) xp, (b) pj^ 
((GeV/c)^), and (c) cos9 for pBe interactions at 530 GeV/c. The cos9 distributions are shown 
separately for each toroid magnet polarity. The integrals of the distributions in (a) and (b) are 
normalized to unity. The integrals of the cos9 distributions within the -0.6< cos6 < 0.6 range are 
normalized to 0.6. The solid and dashed curves in (a) represent empirical fits to the data using 
Eqs. (4.1) and (4.2), respectively. The solid curves in (b) and (c) represent fits using the functions 
shown in Eqs. (4.4) and (4.5), respectively. 
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FIG. 6. Differential distributions for J/tp production as functions of: (a) xp, (b) pj^ 
((GeV/c)^), and (c) cos9 for pBe interactions at 800 GeV/c. The cos9 distributions are shown 
separately for each toroid magnet polarity. The integrals of the distributions in (a) and (b) are 
normalized to unity. The integrals of the cos6 distributions within the -0.6< cosO <0.6 range are 
normalized to 0.6. The solid and dashed curves in (a) represent empirical fits to the data using 
Eqs. (4.1) and (4.2), respectively. The solid curves in (b) and (c) represent fits using the functions 
shown in Eqs. (4.4) and (4.5), respectively. 



29 



T3 

-D 1 



10 



7iBeat515GeV/c 




0.2 



0.4 



0.6 



0.8 



(a) 




10 



-2 



w 
o 
o 



^0.5 



(b) 



^-*± 



't-iq^ 



^^ 



V-^+^^, 



'W 



8 

Pt2 ((GeV/c) 



10 

2x 



(c) 




fait* 




-0.5 



0.5 



cos9 



FIG. 7. Differential distributions for J/tp production as functions of: (a) xp, (b) pj^ 
((GeV/c)^), and (c) cos6 for ir^Be interactions at 515 GeV/c. The integrals of the distribu- 
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FIG. 9. Differential distribution for J/ip production in 800 GeV/c pBe interactions as a 
function of xp compared with: (a) the E771 (pSi) [10] and E789 (pAu) [9] results (each data 
sample is normalized by an integral of the distribution over the 0.0 < xp < 0.135 range); (b) the 
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